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Abstract 

Domain decomposition methods are used for approximate solving boundary 
problems for partial differential equations on parallel computing systems. Spe- 
cific features of unsteady problems are taken into account in the most complete 
way in iteration-free schemes of domain decomposition. Regionally-additive 
schemes are based on different classes of splitting schemes. In this paper we 
highlight a class of domain decomposition schemes which is based on the par- 
tition of the initial domain into subdomains with common boundary nodes. 
Using the partition of unit we have constructed and studied unconditionally 
stable schemes of domain decomposition based on two-component splitting: the 
problem within subdomain and the problem at their boundaries. As an example 
there is considered the Cauchy problem for evolutionary equations of first and 
second order with non-negative self-adjoint operator in a finite Hilbert space. 
The theoretical consideration is supplemented with numerical solving a model 
problem for the two-dimensional parabolic equation. 
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1. Introduction 

Theory and practice of iterative solving stationary boundary value prob- 
lems for partial differential equations is presented comprehensively in the books 
0; HI, lU 3 ■ Different versions of the domain decomposition method with and 
without overlapping of subdomains are used. The approximate solution of un- 
steady problems can be derived via the standard implicit approximations in time 
and solving the corresponding grid problems at the new time level using one or 
another variants of the domain decomposition method for stationary problems. 
Taking into account the transient character of unsteady problems (see, for ex- 
ample, the implementation on the basis of the Schwartz method @, @|), we 
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can construct the optimal iterative methods of domain decomposition where 
the number of iterations is independent of the discretization steps in time and 
space. 

Specific features of unsteady problems are taken into account in the most 
complete way in iteration-free schemes of domain decomposition. In some cases 
it is possible 0, H| to employ only one iteration of the Schwarz alternating 
method for the second order parabolic equation without loss of accuracy of 
the approximate solution. Iteration-free domain decomposition schemes are 
associated with certain variants of the additive (splitting) schemes — regionally- 
additive schemes @- 

Domain decomposition schemes for solving unsteady problems can be classi- 
fied by the method of domain decomposition, choice of decomposition operators 
(exchange of boundary conditions) and used splitting scheme. For differential 
problems it is natural to select domain decomposition methods 

v 

n=|jn Q , n a = n a udn a , a = i,2,...,p (i) 

with overlapping of subdomains (f2 Q| g = U, a n ^ 0) and without overlap- 
ping (Cl a p = 0) [1, HI]. Methods without overlapping of the subdomains are 
associated with an explicit formulation of boundary conditions at the interface 
boundaries. These methods are in common use for solving problems where in 
each particular subdomain its own specific computational grid (triangulation) 
is introduced. To construct homogeneous computational algorithms, domain 
decomposition schemes with the overlapping of subdomains are employed. At 
the minimal overlapping where the width of overlapping is equal to the grid step 
(f2 Q( 3 = 0(h)), domain decomposition methods with the overlapping of subdo- 
mains can often be interpreted as methods without the overlapping of subdo- 
mains supplemented with appropriate boundary conditions of the exchange. 

Domain decomposition (TTJ is associated with an appropriate additive repre- 
sentation of the problem operator: 

p 

A = J2^- (2) 

a=l 

In this case, the operator A a is associated with the solution of some problem in 
the subdomains Q a , a = 1,2,..., p. The most common approach to construct 
the operators of decomposition for solving boundary value problems for partial 
differential equations is based on using of the partition of unit for the compu- 
tational domain. For decomposition (p} we can each separate subdomain Q a 
associates with the function jy a (x), a = 1,2, ...,p such that 

»7«(x) = { ^ 11 %Z «= 1 ' 2 .-^ ( 3 ) 

and also 

p 

^r? a (x) = l, xe!l. (4) 

a=l 
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Suppose, for example, that the operator A is the diffusion operator: 



A = — div fc(x) grad, x 6 SI. 



(5) 



Then for the operators of decomposition we can define the following three basic 
forms: 

Aa = Va A, (6) 

A a = - div k(x)ri a (x) grad, (7) 

A a = Ar) a , a = 1,2,..., p. (8) 

This technique is used beginning with the work [l(| (decomposition ([7J), [H| 
(decomposition ©-(H])), the results of more recent works are summarized in 
the books 0, [l2| ■ Various versions of the decomposition operators correspond 
to using different exchange boundary conditions and ensure the convergence 
of approximate solution in different spaces of grid functions. Special attention 
should be given to issues of constructing decomposition operators for unsteady 
problems with non-selfadjoint operators [lU [3, [lj| . 

For unsteady problems with splitting ([2]) different splitting schemes are used. 



In the theory of additive operator-difference schemes [12, 1|| 17, 18| we need 
to distinguish the case of the simplest two-component splitting. In this case, 
we construct unconditionally stable factorized splitting schemes, such as the 
classical scheme of alternating directions, predictor-corrector scheme. Two- 
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component regionally-additive schemes are constructed and studied in 
[l|| as well as in the above papers IH 3 HI for convection-diffusion problems. 

In application of domain decomposition methods the splitting of problem 
operator into the sum of three or more non-commutative operators (p > 2 in 
([2])) is of great interest. Classic schemes 2 12, 18 1 of multi-component splitting 
are based on the concept of summarized approximation. Additively-averaged 
schemes of summarized approximation 



are more explicitly oriented to 
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parallel computations. Regionally-additive schemes of component-wise splitting 
are investigated in 21]. A variant of two-component splitting with the Crank- 



Nicolson scheme for the individual subproblems with the minimal overlapping 
and decomposition ([7]) is considered in the article [22 1. 

Nowadays, the schemes of full approximation are in common use for the 
general multi-component splitting. In this regard, we note regularized additive 
schemes [23| where the condition of stability is achieved due to perturbations 
of operators of the difference scheme. In the vector additive schemes 24| 25| 
instead of one equation we solve a system of similar equations. Such schemes are 
also constructed for the evolutionary equations of second order 2(| 27l. V ector 
regionally-additive schemes are investigated in 28, 29|. In the work [3d] there 
are proposed more general regularized schemes of domain decomposition with 
different structures for both the splitting operators and operators of the grid 
problem at a new time-level. 

Among other domain decomposition methods for solving boundary value 
problems for parabolic equations it is necessary to hi ghl i ght ex p licit- im plic it 
methods considered in many papers (see, for example, 31, 32, 33[ 34, HI, 36|). 
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The domain decomposition in this case is performed without the overlapping of 
subdomains and the transition to a new time-level is organized as follows. First, 
the approximate solution at the common boundaries of subdomains is predicted 
using the explicit scheme. Next, these boundary conditions are used to derive 
the approximate solution within individual subdomains. And finally, a correc- 
tion of interface boundary conditions is carried out using implicit schemes. It 
will be shown below, that such schemes of domain decomposition are completely 
fit in the above general scheme of decomposition methods at a special domain 
decomposition with the choice of operators according to decomposition ([5]) . 

In this paper we construct domain decomposition schemes for parabolic and 
hyperbolic equations with self- adjoint elliptic operators of second order. Un- 
conditionally stable factorized regionally-additive schemes are constructed using 
decomposition (J6j> and the two-component and general multi-component split- 
ting. Domain decomposition schemes with a self-adjoint operator for the grid 
problem at the new time-level are derived, that allows to construct on their 
basis iterative methods of domain decomposition-conjugate gradients for solv- 
ing stationary problems. The paper is organized as follows. In Section 2 we 
formulate the model Cauchy problem for parabolic and hyperbolic equations 
in a rectangle. Next, Section 3 provides the stability conditions for the stan- 
dard two- and tree-level implicit schemes with weights for model problems in 
a rectangle. The domain decomposition and construction of the operators are 
discussed in Section 4. The possibilities of the standard factorized schemes for 
domain decomposition are considered in Section 5. The stability condition, ap- 
propriate a priori estimates for the approximate solution and convergence rate 
estimate are derived for factorized regionally-additive schemes. In Section 6 
there are constructed multi-component splitting schemes. Hyperbolic equations 
of second order are considered in Section 7. Theoretical results are illustrated 
by the numerical results presented in Section 8. 

2. Model boundary problems 

Let us consider a model boundary value problem for the parabolic equation 
of second order. In a bounded domain Q the unknown function it(x, t) satisfies 
the following equation 

lr fab (•«■■>£) -/(-■<>. "<'<- T - < 9 > 

a — 1 x ' 

where fc(x) > n > 0, x G ft. Equation © is supplemented with the homoge- 
neous Dirichlet boundary conditions 

u(x, t) = 0, x e <9ft, < t < T. (10) 

In addition, the initial condition is prescribed 

w(x,0) =«°(x), x6!l. (11) 
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Unsteady diffusion problem (|5|)- (fTTj) is considered on the set of functions 
u(x,t) satisfying boundary conditions flTUl) . Then instead of ([9]), (fTUl) we use 
the operator-differential equation 

^+Au = f(t), 0<t<T. (12) 



The Cauchy problem is considered for evolutionary equation (|12[) : 

u(0) = u°. (13) 
For the diffusion operator we assume 



^=-£^( fc(x) £) 

a— 1 v ' 



On the set of functions (fTO)) let us define the Hilbert space % = £2^) with 
the scalar product and norm 

(u, v) = / it(x)v(x)<ix, ||u|| = (it, u) 1 / 2 . 



In H the operator of the diffusive transport A is self-adjoint and positive definite: 

A = A*>k6E, 6 = 5(Q)>0, (14) 

where £ is the identity operator in H. 

We present now the simplest a priori estimate for the solution of problem 
(fT2")) - (fT4")) which will be for us the check point for the considering grid problems. 
The self-adjoint positive definite operator V can be associated with the Hilbert 
space T-Lv having the inner product and norm 

1 /2 

(u, v)v = (Vu, v), \\u\\t> = (u, u) v 

respectively. In H multiply scalar ly equation (fT2"j) by Au. In view of p^|) we 
obtain inequality 

' d '\u\\ 2 A +\\Au\\ 2 = (f,Au). (15) 



Taking into account 
from (|15|) we have 



2 dt 

(f,Au)<\\Au\\ 2 + \\\ff, 



i t M\<\\\ff. 
In view of the Gronwall lemma we obtain the desired estimate 

H3i<ll« ll3i+ fwmfde, (i6) 
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which expresses the stability of the solution of problem (fT2)) -(fIl | with respect 
to the initial data and right-hand side. 

In addition to the parabolic equation ©, we consider the hyperbolic equa- 
tion 



-'«"■'>• ««• 



(17) 



with boundary conditions (1101) . Equation (|17p is supplemented with two initial 
conditions 

u(x,0)=u°(x), ^(x,0)=<; (x), xeO. (18) 

Problem pH|) . (fTTjl . (|T5)) is associated with the following Cauchy problem for 
the evolutionary equation of second order: 



— +Au = f(t), 0<t<T, 

u(0)=u°, ^{0)=v°. 
at 

Multiply scalar ly equation (TK)|) by Adu/dt and obtain 

2 \ 



ld_ 

2 dt 



du 
~dt 



\Au\ 



(19) 
(20) 



For the right-hand side we use the estimate 



, du 
dt 



< 



du 
~d~t 



+ ^ll/llV 



The result is 



where 



f t \\u\\l<\\u\\t + \\f\\% 



mi; 



du 
~dt 



\Au\ 



The desired a priori estimate 



\\u{t)\\l < exp(t) \\AuY + WvYa + eM-O)\\f(0)\\AdO (21) 



expresses the stability with respect to the initial data and right-hand side of the 
Cauchy problem for operator-differential equation (jlT 
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3. Standard difference approximations 

We will conduct a detailed study of approximations in space and time using 
as an example the boundary problems in a rectangle 

fi = { x | x = (xi,x 2 ), < x a < l a , a = 1, 2}. 

The approximate solution is given at the nodes of a uniform rectangular grid Q: 

U) = {k\x = (xi,x 2 ), x a =i a h a , i a = 0,l,...,N a , N a h a = l a } 

and let oj be the set of internal nodes (Q = uj U dm). For the grid functions 
y(x) = 0, x 6 du we define the Hilbcrt space H = L 2 (co) with the scalar 
product and norm 



{V,w) = ^2y(x)w(x.)h 1 h 2 , \\y\\ = {y,yf 



/2 



Assuming that the coefficient fc(x) in SI is sufficiently smooth, we take the 
grid operator of the diffusion as 



1 

1 

hi 
1 

h. 



Ay = - jok(xi + 0.5/ii, x 2 )(y(xi + hi,x 2 ) - y(xi,x 2 )) 



+-^k(xi - 0.5h i ,x 2 )(y(x 1 ,x 2 ) - y(xi - hi,x 2 )) 



2 k{x\,x 2 + 0.5h 2 )(y(xi,x 2 + h 2 ) - y(xi,x 2 )) 



+ j^k(xi,X2 - Q.bh 2 )(y{xi 1 x 2 ) - y(xi,x 2 - h 2 )). (22) 
In H the operator A is self-adjoint and positive definite: 

A = A* > k(6i + S 2 )E, 5 a = ^sin 2 ^, a = 1,2. (23) 

hi 2l a 

After the approximation in space we go from ([9]), (fTO]) to the differential- 
difference equation 

^+A V = /(x,t), xe^, 0<t<T. (24) 

Taking into account (fTTj) . let us supplement equation (|24|) with the initial con- 
dition 

y(x,0)=tt°(x)» xe^. (25) 



For the solution of the differential-difference Cauchy problem (|24|) . ([25]) the 
following a priori estimate holds (see (|16[) ) 

IMIi<U«°lli + twmfdo- (26) 

Jo 
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Similarly, the approximation in space leads us from (fTU|) . (|T7|) . (|T5|) to the 
problem 

+ Ay = /(x,t), xew, 0<i<T, (27) 

l/(x,Q)=« (x) ) ^(x,0)=«°(x), xew. (28) 
The grid analog of (f2~Tj) is the estimate 

||, /( /)||- .x,,(/) I |L4 U °|| 2 +II«°IIa + / exp(-0)[[/(0)||^0 I . (2«)i 



where 

2 



lyll* 



A 



The emphasis now is on the approximation in time. In the construction of 
domain decomposition schemes for problem (|24|) . (|25|) . the starting point for 
us is the usual two-level schemes. Let r be of a uniform time-step and let 
y n = y(t n ), t n = nr, n = 0, 1, ...,N, Nt = T. Equation ((24]) is approximated 
by a two-level scheme with weights 

* ^ + A( ( ry"+ 1 + (l- ( 7) 2 ;») = ^, n = 0, 1, .., N - 1, (30) 

T 

where, for example, <p n = f(at n+1 + (1 —a)t n ). It is supplemented by the initial 
condition 

y° = u°. (31) 

Difference scheme (j30|) , ([3Tj) has the approximation error 0(t 2 + — 1/2)t+/i 2 ), 
where ft 2 = (fr 2 + h 2 )/2. 

Theorem 1. Difference scheme i30\) , h31}) is unconditionally stable for a > 
1/2, and for the numerical solution the estimate 

\\y n+1 \\ 2 D<\\y n \\l + ^y n \\ 2 , n = 0,l,...,iV-l, (32) 



holds, where 



D= .l-i- i a - i ) r_l ;J 



Proof. Let write difference scheme (1301) as 



and multiply scalarly it by TA(y n+1 + y"). Using the fact that a > 1/2 the 
operator D > A, we have 

II2/" +1 IId - h n \\l + \U{y n+l + y n )f = <v\ A( y n+1 + y n )). 
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Taking into account 



(p* A(y n+1 + y n )) < \\\A{y n+1 + y n )\\ 2 + i||^"|| 2 , 



we obtain the required estimate ([32 



A priori estimate (|52"j) for the solution of problem (j3"0)) . (I5T|) is a grid analog 
of the a priori estimate (|26[) for the solution of differential-difference problem 

d2!, (GH> Cd = a + o(t)). 

To solve numerically problem (|27|) . (f28]h it is natural to use three- level 
schemes of second order accuracy in time. Let 



n i 1 



+ A(aj/ n+1 + (1 - 2a)y n + ay^ 1 ) = <p n , 



n=l,2,...,N-l, 



(33) 



where, for example, ip n — f(t n ). In view of ([25)1 we can for the solution of 
equation (|27|i approximate the initial condition as follows: 



o o y 1 -y" 

V = u , 



= v° + -Au°). 

T 2 



(34) 



The error of difference scheme (|33"]l. (|34]) is 0(t 2 + h 2 ). 

Theorem 2. Difference scheme A33\) , is unconditionally stable for a > 

1/4, and for the numerical solution the estimate 



exp(r) 



2 exp(0.5r) - 1 



<p n \\\, n = 0,l,...,N-l, (35) 



holds, where 



S n = 



y n _ y n-l 



A 



y n + y n-l 



D = A 



r 2 A 2 . 



Proof. We introduce the notation 



n + n-1 

C = — — — , v 



y n _ y n-l 



Taking into account the identity 



y n = \{y n+1 + 2y n + y^ 1 ) - \{v n+1 - 2y n + y^ 1 ), 



ay n+ ' + (1 - 2a)y n + ay n ~" = y n + a(y n +' - 2y n + y^') 
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we rewrite (|33p as 

(*+(,_ J) ^)*^ + ji <^t£_^. (36) 

Multiply scalarly (|36|) in H by 

2A(C™ +1 - C") = TA(n n+1 + rf). 
With this notation for a > 1/4 we obtain 

S n+1 - S n = T(<p n ,A(ri n+1 +r) n )). (37) 
Using the estimates for the right-hand side 

TA(<p», ( V n+1 + r, n )) < y £ \\fj n+1 + V n f A + T -e\\ V n \\\, 

||r/" +1 +ry"||i<2(||ry"+ 1 ||i + ||ry"||i), 
with e > 0, from (l37l) we obtain 

(l_I)^+i<(l + l) 5 " + ^||^||i. (38) 

We choose e so that 

T 

1 = exp(— 0.5t), 

and therefore 

1 + T - = exp(0.5r). 

With this in mind from (f3"5)) we obtain the level- wise stability estimate ([35]) . 

Estimate (1331) can be treated as the grid analog of the a priori estimate (|2U)) . 
For difference schemes (j3"0)) . ([3T|) and ([33]). ([3^]) we can obtain many other a 
priori estimates of stability with respect to the initial data and right-hand side 
[Si, [l|| . We have restricted to only those estimates that we can associate with 
the corresponding estimates for the domain decomposition schemes considered 
below. 

4. Substructuring domain decomposition 

Let us consider a special class of domain decomposition methods. At the 
discrete level we define a set of interface nodes inside the domain and then solve 
the subproblems separately inside the subdomains. At the continuous level, this 
decomposition is associated with subdomains the width of which is equal to the 
corresponding discretization step in space. We illustrate our consideration on 
the model grid problems in a rectangular. 

The computational grid ui is partitioned into rectangular subdomains of a 
coarse grid with the step h. The boundaries of the subdomains (direct lines) 
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Figure 1: Grid decomposition 



consist of the nodes of the fine computational grid. Denote this set of inte- 
rior boundary nodes as Q. A fragment of the grid is shown in Fig. [T] Such 
a decomposition of the fine computational grid corresponds to the domain de- 
composition depicted in Fig. [2] fl = fix U H2, Q12 = 0- Subdomain VL2 is a 
lattice, width of the individual edges of the lattice is h. Domain 51i consists of 
disconnected individual subdomains. 

The partition of unit for ([3]), Q we associate with the corresponding additive 
representation of the identity operator E in the space of grid functions H, 
defined on the set of internal nodes of w. Let 

p 

Y^Xc^E, Xa >0, a = 1,2,..., p. (39) 

Q = l 

Similarly ©, the operators of decomposition can be given in the form 

A a = X aA a = 1,2,..., p. (40) 

In view of (|39p. in this splitting we have for the problem operator the following 
additive representation 

A = J2A a . (41) 

a=l 
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Figure 2: Domain decomposition 



Splitting (|4Tj) allows us to go from equation ([24]) to the equation 



dy 
dt 



^2A a y = f(x,t), xew, 0<i<T. 



(42) 



a = l 



Direct construction of various splitting schemes for problem (l25l , (|42j) is com- 
plicated by the fact that individual operator terms A a , a = 1, 2, ...,p do not 
inherit the basic properties of the operator A — the self-adjointness and non- 
negativity. However, using decomposition operators (14*(J1) . equation (|4*21 can 
be easy transformed in the symmetric form. Multiplying equation (1421) by the 
self-adjoint operator A, we obtain the equation 



a— 1 



(43) 



where the operators 

B = A, A a = AxaA, a = l,2,...,p 
are self-adjoint and non-negative. Moreover, we can introduce new variables 
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v = A 1 / 2 !) and instead of f4"3")l we can consider the equation 

C ^ + ^A a v = A 1 l 2 f{x,t), xew, 0<i<T, (44) 

a=l 

with self-adjoint and non-negative operators 

A a = A 1 ' 2 Xa A 1 ' 2 , a = 1,2,..., p. 

Standard estimates for the solution of equation (35)1 in the norm of H (for \\v\\) 
correspond to using estimates in Ha (for ||j/||.a)< This explains our unusual in 
some sense choice of the the priori estimate (j26|) for problem (p4|) . (|25|) and 
estimate (25} for problem ([27]). (25]). 

The particular specification of the decomposition operators of type (|39[) . 
P0|) is provided via the selection of terms x Q , a — 1?2, ■ •-,£>■ Some advanced 
features are discussed below, but we start from the simplest version. If we use 
substructuring domain decomposition (see Fig. [T]), it is natural to put 

X2(x) = j^ Xl (x) = l- Xa (x), xec. (45) 

The operator A2 is associated with interface nodes Q, whereas A\ — with the 
internal nodes of subdomains. 



5. Factorized schemes of domain decomposition 

After selecting the operators in decomposition (|4Tj) the construction of do- 
main decomposition schemes is carried out using one or another additive schemes. 
For (|40)) . (|45|) we can consider the simplest two-component (p = 2) splitting 
schemes. In this situation, we can try to use the operator analogues of the 
classical schemes of alternating directions H, 39[ . 

We begin with the scheme of stabilizing correction [3^] , where the transition 
to a new time level in problem (|23|) . (|4"2")l with p — 2 is performed as follows: 

v n+l/2 _ n 

V - U — + A lV n+1 / 2 + A 2 y n = <p n , (46) 

T 

n+1 _ n 

U - U — + A lV n+1 / 2 + A 2 y n+1 = p n , (47) 

T 

where, for example, ip n = f(t n+1 ), n = 0, 1, ...,N - 1. Equations (46]), are 
complemented by the initial condition pip . 
If decomposition (|4T)t , (gSJ) is used, we have 

„,n+l/2 _ n 

U - y — + X iAy n+1/2 + X 2Ay n = <p n , (48) 

T 

V - V - + X iAy n+1/2 + X 2Ay n+1 = V n . (49) 
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The implementation of this scheme can be different. 
Taking into account that 

<f n = Xi<P n + X2<p n , 

let us introduce the auxiliary function y n+1 / 2 and divide equation f)48[) into two 
ones: 

~n+l/2 _ n 

y - »— + X 2Ay n = X2 l fi n , (50) 

T 

,,"+1/2 _ ^n+1/2 

y - v - + xiAy" +1/2 = xw n - (5i) 

r 

The function y™+ 1 / 2 is determined via the explicit scheme (|50l) . Moreover, 
taking into account (j45|). the calculations are performed only on the set of 
interface nodes. 

Stage 1. Evaluation of the conditions at the boundaries of the subdomains 
via the explicit scheme: 

~n+l/2 _ n 

y — + Ay n =y n , xe£, 

T 

r +1/2 = y n , x^s. 

After such a predictor of boundary conditions we solve problems in subdomains 
GB- 

Stage 2. Evaluation of the solution in subdomains using the implicit scheme: 

y v + Ay n+i/2 = » x g ~ 

T 

y»+V2 = yn+V2 5 x g £. 

The last step is to correct conditions at the boundaries, which provides, in 
particular, the stability of the approximate solution. For the subdomains it is 
convenient to replace equation (|4"5j) by the difference of (|4"5|) . (|4"gj) : 

, n+l _ n+1/2 

* 2 + X 2A(y n+1 - y n ) = 0. 

T 

Taking into account (|45p , we calculate the approximate solution at the new time 
level. 

Stage 3. Correction of the conditions at the boundaries of subdomains via 
the implicit scheme: 

n+l _ n 
T 

y n+1 = y n+1/2 , x^S. 

This numerical implementation (stages 1-3) of regionally-additive scheme 



(|4"5|) , (|48p , P^)l is nothing but the scheme of the domain decomposition 3l|, [32 
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IH 34, 2,5, 36]) with the explicit- implicit procedure for calculating the boundary 
conditions at the boundaries of subdomains. 

Regionally- additive scheme (|45p . ([4"8")) . f4"5|) has the first order approximation 
in t. It is possible to use the schemes of second order where 

H+l/2 _ n 

V - — JL + X iAy n+1 ' 2 + X2 Ay n = v », (52) 

t/2 

n+l _ n+1/2 

y - f- + xiAy n +^ 2 + X2Ay n+1 ^ip n . (53) 

t/2 

with <p n = f(t n+1 / 2 ). Schemes ggj, 0SJ and ([52]), (J53J we consider as the 
operator analogs of the classical schemes of alternating directions. They are 
special cases of more general factorized schemes. 
Consider the factorized scheme 

n+l _ n 

B X B 2 » U - + Ay n = <p n , (54) 

T 

where 

B a = E + <TT Xa A, a = 1,2, (55) 

with the right-hand side specified in the form ip n = f((rt n+1 + (1 — a)t n ). Direct 
substitutions verify that scheme (|5i|) . (|55l) coincides with scheme (|4*51) . P^l) at 
ex = 1 and with scheme ([521 . ([531) at cr = 1/2. 

For the factorized scheme ([541 . ([55)) it is possible to use the three-stage 
computational implementation with explicit-implicit calculations of interface 
boundary conditions. We introduce, for example, the new grid function y n+1 
and instead of ([54| in view of ([55]) we solve two differential equations: 



~n+l _ n 

{E + cjtXiA) U - U — +A y n = <p n , (56) 

T 

n+l _ n rrn+1 _ n 

(E + arxiA) V - V - = V - V - . (57) 

T T 

Taking into account ([45[) . we obtain from (|56|) that for nodes at common bound- 
aries {Stage 1 — the explicit scheme for boundary nodes): 

yn+l _ n 

y — + Ay n = ip n , x e Q. (58) 

T 

For subdomains we have: 

n+l _ n 

{E + arA)^ y — + Ay n = tp n , x 4 Q. 

T 

This corresponds to (Stage 2 — the implicit scheme in the subdomains) usage of 
the implicit scheme with weight a for the difference solution in the subdomains. 
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The implementation of (|57|l (Stage 3 — the implicit scheme for the boundary 
nodes) in view of (|45|) is: 

(E + arA)- — = — , xe£, (59) 

T T 

y n+1 = r +1 , x^S. 

Taking into account (|58|) . equation (f59| can be written as 

n+l _ n 

(E + arA)^ y —+Ay n = ip n , xew. 

T 

At this stage all computational work is associated only with correction of the 
internal boundary conditions via this implicit scheme with weights. 

Theorem 3. Factorized regionally- additive difference scheme LS9\) , \5$ , \55)) 
is unconditionally stable for a > 1/2, and for the difference solution the follow- 
ing estimate holds 

\\B2y n+1 \\ A < ||S 2 y n |U+r||BrVlU, n = 0,1, N — 1. (60) 



Proof. It is convenient firstly to symmetrize the factorized scheme (I54|) . ([55]) . 

Let v n = A 1 / 2 y n and 

B a = E + arA a , A a = A 1 l 2 Xa A 1 l 2 , a = 1,2. 
Then equation (|54[) can be rewritten as 

B X B 2 + Av n = A 1 / V- (61) 

T 

Assuming that Biv n = io n , from (|61l) we obtain 

w "+l = Sw n + tB^A 1 ' 2 ^ 1 , (62) 
where the operator of the transition to the new time level 

S = E-tB^ABz 1 . (63) 
Taking into account the above notation, from (|S3l we obtain 

S=^-±E+ 1 ^B^ 1 (B 1 B 2 -2a(A 1 +A 2 ))B^ 2 
la la 

= — ^ E + S1S2, 

la 

where 

S a = (E + aTA a )- 1 (E-aTA a ), a = 1,2. 
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If a > 0, taking into account the non-negativity of the operators A Q , a = 1, 2, 
we have 

||S a ||<l, a =1,2. 

With stronger restrictions a > 1/2 we find that ||5|| < 1. From (|52p we obtain 
the estimate 

= IKIl + TiiBr^Vll- 

This is the required estimate 

The fundamental issue in the construction of domain decomposition schemes 
for unsteady problems is to estimate the convergence rate for the approximate 
solution. Accuracy depends on a computational grid (the width of the overlap- 
ping) and therefore regionally- additive schemes belong to the class of condition- 
ally convergent. The situation can be illustrated by the example of the above 
factorized decomposition schemes (|54l) . (1551) . 

Analysis of the accuracy will be conducted in the standard way by consid- 
ering the corresponding problem for the error 

z n (x)=y"(x)-u«(x), xGc, 

where M n (x) = u(x, t n ) is the exact solution of the differential problem (|9T)~(fTTj). 
From (|54l) . (1551) we obtain the problem for the error 

7 n+l n 

Si B 2 - + Az n = V/ 1 , (64) 

T 

z° = 0. (65) 
In view of ([60]) for problem (|64|) . (|65|) we have 



\\B 2 z n+l \\ A <Y,r\\B^^ k \\ A , n = 0,l,...,N-l. (66) 

fc=0 

For the approximation error we have 

? .n+l _ u n 

^"=^"-5x^2 Au n . (67) 

r 

Taking into account (|55|) . from (|67|) we obtain 

„ ( / 1\ A M n+1 -u n A 
= <p n - I E + I a - - J rAj A- 

^ = -a 2 T 2 Xl A X 2A 
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The first term of the error is the standard one for the schemes with weights, 
whereas the second term results from the splitting of subdomains. For suffi- 
ciently smooth solutions of problem (|§))- (fTTj) we have 



Let us consider the term tp^ m more detail. 

Taking into account (|66l) and introduced in the proof of Theorem [3] notation, 
we have 

HBrV?iu = iiflr 1 ^ 1/ V n ii 



^2 2 

a t 



B^A^xxAx-zA 



u n+l _ u n 



,n+l 



< err 



A^A- 



where 
Thus 



Q = {E + arA 1 y 1 arA 1 . 



\\B^\\ a = 0(<tt\\ X 2\\a). 
These arguments allow us to formulate the following statement. 

Theorem 4. For the error of the factorized regionally- additive difference scheme 
H39\) , \5J$ , A55\) with a > 1/2 we have for problem Hfy- llll]) the following esti- 
mate 

\\B 2 z n+1 \\ A <M^h 2 + T 2 + ^-^r + arllxallA (68) 

For considered here substructuring domain decomposition schemes with the 
grid elliptic operators of second order (f2"2"]l and splitting , (|4"5]) estimate (j6"5j) 
gives 

\\B 2 z n+1 \\ A <M(h 2 +r 2 + L-^jr + arh-^h- 1 / 2 ] . (69) 

Note also that the use of the scheme with a = 1/2 does not increase the order 
of accuracy. But in this case the main error term is two times lower compared 
to a = 1. 

A slightly different algorithm can be implemented. Instead of we apply 
the factorized scheme 



n+l _ n 

B 2 B X U - y -+Ay n = <p n , 



(70) 



i.e. we comutate operators B\ and B 2 . 

The implementation of scheme (I7D1 will differ slightly from the implementa- 
tion of scheme (l54l). Similar to (1561). (15711 we have 



(71) 



(E + utxiAY- y — + Ay n = <p« 
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(E + arxiA) * U — = U - V — . (72) 

T T 

At stage (|7T|) we use the implicit scheme for the nodes at the boundaries of the 
subdomains and explicit scheme in the subdomains. Note that for the explicit 
scheme it is enough to evaluate only the boundary nodes. At stage (|72|) the 
solution in the subdomains is calculated using the implicit scheme. Thus, the 
computational cost in case of the factorized scheme ([70)1 remains practically the 
same as for scheme (l54l). 



6. Schemes of multi-component splitting 

Constructed above factorized schemes of the two-component splitting can be 
generalized in various directions. The most fundamental issue is to construct 
such schemes in the case of general multi-component splitting. 
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Figure 3: Three-component decomposition without the overlapping of subdomains 



The need for such an extension results from, in particular, calculations of 
conditions at the boundaries of subdomains, i.e. the solution of problems on 
graphs for two-dimensional problems. In the considered two-dimensional prob- 
lems in a rectangle and rectangular grids, the implementation of, for example, 
P5|) does not face significant problems. However, for more general situations, 
for example, for three-dimensional boundary value problems, the solution of 
these grid problems can be difficult. Such considerations lead us to the need 
of constructing procedures of decomposition for the set of boundary nodes of 
subdomains. A characteristic example is shown in Fig. [3] The set of boundary 
nodes is divided into two parts: lu = uj s U u} m . Here the set of nodes at the 
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Figure 4: Three-component decomposition with the overlapping of subdomains 



boundary of two subdomains is denoted as u s (in Fig. [3] it is depicted as •). 
The set of nodes that lie at the boundaries of a greater number of subdomains 
is designated as u m (in Fig. [3] it is presented as □). 

Instead of the two-component splitting (|39|) . (j45|) . we use now the three- 
component splitting (|39|) with p = 3 and 



, , J 1, X e U s , , , J 1, X e U) m , 

X2(X) = 1 0, x*Q„ * 3(X H 0, x*£5 m> 

Xi(x) = l-% 2 (x) -x 3 (x), xew. (73) 

With such a decomposition calculations in some parts of the sub-domain bound- 
aries (on the set u s ) can be performed independently using known conditions 
at the nodes of crossing (on the set u m ). 

Local computations of of the solution at boundary crossings introduces ad- 
ditional errors. To improve the accuracy of the approximate solution at the 
boundaries of subdomains, it is possible to apply algorithms with the overlap- 
ping of subdomains. Such a situation at the grid level is shown in Fig.|4] There is 
highlighted the set of boundary nodes u m , which lie near the boundary crossing 
and u s n u rn ^ 0. With this in mind, instead of (|73[) we set 

, , / > 0, x G u s , , , / > 0, x 6 u m , 

X2(X) = \0, x*D„ * 3(x Ho, x*C m> 

X2(x) +% 3 (x) = 1, xew, xi(x) = 1 - X2(x) - X3(x), xew. (74) 
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For the general multi-component (p > 2) decomposition it is possible to 
construct in a more simple way regularized additive schemes [l2l . |23| . For solving 
problem ([23]) . (HHJ), (H2"]l we can use the additive scheme of full approximation 

n+l _ n 

U - U - + Ay n = <p n , (75) 

T 

where 

v 

A = J2A a , A a = (E + <jT Xa A)- 1 x a A 1 a = 1,2,..., p. (76) 

a=l 

This scheme is characterized by the fact that each operator term XaA, a = 
1, 2, is perturbed with an error 0{t). 

Theorem 5. Regularized difference scheme J75| ), J76'| ) is unconditionally stable 
for a > p/2, and for the difference solution we have the estimate 

\\y n+1 \\A<\\y n \\A + r\\<p n \\ A , » = 0,1,...,JV-1. (77) 

Proof. The operator A can be written in the form 

A = Y j A~ 1 / 2 R a A 1 '\ (78) 

where 

R a = {E + aTC a y 1 C a , C a = A 1 / 2 Xa A 1 / 2 , a = 1,2,..., p. 

Thus 

C«=C*>0, R a =R* a >0, a = 1,2,..., p. 
With this in mind the difference scheme (f75|) . ([7o]l is written as 



J2 R * vn = A 1/2 tp n , (79) 



a=l 

where, as before, v n = A x l 2 y n . From (j?9")) we have 

v n+1 = Sv n + TA 1 > 2 <p n (80) 
with the transition operator 

v 



S = E - tJ2 r *- 



a=l 

Using this representation, we set 

p 



1 

S=~y^S a , S a =E-prR ai a = 1,2,..., p. (81) 
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For the individual terms with the above notation we obtain 

S a = (E + aTC a )- l (E + (a - jj)rC a ). 

Under the constraints er > p/2 we have \\S a \\ _ 1, which allow us to obtain from 
(|5D|) the following estimate 

II^HIKII+rii^/Vll, 

which is nothing but (|77p. 

Standard finite-difference schemes of component-wise splitting lij, 17, [l8| 
can be easy constructed using the transition operator. With regard to our prob- 
lem, we shall again start with notation (T5U|) . but instead the additive structure 
(see (|8Tj) employ multiplicative one 

l 

5=[[5o, S a = E-rR a , a = 1,2,..., p. (82) 

a— p 

In this case we have \\S a \\ < 1 for a > 1/2. 

The implementation of component-wise splitting scheme is performed as a 
sequence of intermediate difference problems similar to (|80p : 

w «+«/p = 1 5 aU "+(«-l)/P +ry il/2 ¥ ,n ) a = 1,2,..., p. (83) 
Comparing with (JSUJ) , (|8ip . we obtain 

1 a-l 

v n = e n ^s- 

a=p /3=p 

Without loss of accuracy, we can consider only the simplest choice for <^™, a = 
1, 2, where 



With this notation the difference equations (1831) can be written as follows 

„ n+a/p _ n+(a-l)/p 

V - V - + (E + aT Xa A)- l Xa Ay n +^/v = 

T 

a = 1,2,..., p. (85) 
We can formulate now the following statement. 

Theorem 6. Additive component-wise splitting scheme fl£^[ ), h85\) is uncondi- 
tionally stable for a > 1/2, and for the difference solution estimate |77|) holds. 
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The considering scheme has the first-order approximation in time. However, 
in the case of two-component splitting at a = 1/2 the approximation error 
is 0(t 2 ) (See, e.g., [ioj). This variant is used in [22| to construct domain 
decomposition schemes. For the general multi-component splitting the additive 
schemes of second order in time are based on the symmetrization of transition 
operator [4l|,|42[. In this case, instead of (1521 we can use, for example, 

p i _ 1 
5= n^Il 5Q ' S «=(E+~C a ) (£?-~C Q ), a -1,2,..., p. 

/3=1 a=p 

Thus we make two half-steps in time for a > 1/2 in sequential solving problems 
for operators XaA, a — 1; 2, ...,p, and then for the operators X/3^j = P,P — 
1,...,1. 

Regularized scheme (|75|) . (|76|) can be written in the form similar to (|85|) : 



,i+cx/p _ n+(a-\)/p 



a = 1,2,..., p. (86) 

In contrast to (|55|) here the obtained y n +( Q ~ 1 )/p is used for solving the problem 
for y n + a /P only partially. This increasing of explicitness results in a more strong 
condition of stability (instead of a > 1/2 we have a > p/2). 

The numerical implementation of scheme (1841) , (|85p is shown below for de- 
composition P5|) . In accordance with (|8"4")l we set 

a = 1,2,..., p. 

From ([85]) we obtain 

„.n+l/2 _ n 

(E + ar X iA) U - u — + Xl Ay n = (E + ar X iA) XaV n , (87) 

r 

n+l _ ri+1/2 

(£ + VTX2AY- U - + X 2Ay n+1 ' 2 = (E + ar X iA)x^ n . (88) 

T 

In finding y"+ 1 / 2 from (|87[) (Stage 1) we solve boundary value problems in 
subdomains using the implicit scheme. The boundary conditions are taken from 
the previous time level, i.e., 

n+l/2 _ n 

(E + gtA)^- y — + Ay n = (E + <rrA)<p n , x Q. 

T 

^+1/2 =y n x£ ^ 

Conditions at the common boundaries are corrected during evaluation y n+1 
from dFJ) (Stage 2): 

n+l _ n+l/2 

(E + gtA)^- y - + Ay n+1 ' 2 = (E + (jTA)tp n , x e u>, 
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y n+1 = y n+1/2 , x^S. 

The numerical implementation of the component-wise splitting scheme (|84p , 
(|55|) is slightly reduced in compare with the factorized domain decomposition 
scheme (|4*5j) . (|4*5)) (there is no explicit calculations of the interface boundary 
conditions). Similarly, two-stage implementation takes place for regularized 
scheme (|H3D, dHSD- 



7. Hyperbolic equations of second order 

Possibilities of constructing domain decomposition schemes to solve bound- 
ary value problems for hyperbolic equation of second order (ITU)) , (TTT)) . (IT51) arc 
more restricted. Here we note only the regularized schemes similar to (|75[) . (|76l) 
for the parabolic problem (|9"j)-([TT]), 



n i „.n— 1 



(89) 



where 



A = Y^A a , A a = {E + aT 2 Xa A)- 1 Xa A 1 a =1,2,..., p. (90) 

This scheme has the second order of accuracy in time. The following statement 
is true. 



Theorem 7. Regularized difference scheme {3$ , &89\) . i9ty) is unconditionally 
stable for a > p/4, and for the difference solution the following estimate is 
satisfied 

S n+1 < exp(r)5" + y — g*^- - n = 1,2, N — 1, (91) 

where 



y n _ y n-l 


2 


y n + y n-l 


2 


T 


+ 

D 


2 


AA 



D = D* = A \ E A 

4 



Proof. The proof is conducted similarly to Theorem[2] Similarly (l78l) . for the 
operator A, taking into account (|90|) , we have representation 



A = ]T A^^R^A 1 ' 2 , 



(92) 



where now 

R a = {E + oT 2 C a )- l C a , C a = A l ' 2 Xa A 1 ' 2 , a = l,2,...,p 
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with self-adjoint and non-negative operators C a R ai a = 1, 2, 
Difference scheme (l89|) , (|90|) can be written as 



R a v n = A 1 ' 2 



for v n = ^l 1 / 2 !/". Using the notation 

c n = 

write ([M]) as 



E ~ t r )—r L+R —Y^ 



where 



R — R — ^ Ra- 



a=l 

Multiply dMl) scalarly in H by 
and obtain the equality 

S n+1 _ gn _ r ^l/2y,» 5 (^"+1 + f? «)) j 

where 

^-||'7 n || 2 5 + llC"H 2 « 



D = E — —R. 

4 

and D > 0. For the first term we have 



y n _ v n-l 


2 

+ 


v n +y n-l 


T 


D 


2 



D" - D™- 1 


2 


y" - y n ~ x 


2 


T 


D 


T 


D 



that results from 

A 1 / 2 DA 1 / 2 =A- T 1a 1 ' 2 RA 1 ' 2 = aU- ^Aj = L 

To prove the inequality D > for a > p/4, we set 

1 p ~ t 2 
D=~y^D a , D a = E -pR ai a = 1,2,..., p. 

^ a=l 
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For each individual term we have D a , a = 1, 2, ...,f>, if 
£7 - — p(E + oSCa)- X C a > 0. 

We have 

E + ar 2 C a - -jpC a > E. 

for a > p/4 for each a = 1, 2, 

For the right-hand side of ([55)1 with <£>™ = A x / 2 (p n , we use estimates 

r(fc (rj n+1 +rf)) < ^h n+1 +v n \\% + ^ll^lll- 

|| 7? »+i +7? «|||<2(||rr+ 1 ||| ) + || ? ni| ) ). 
We obtain the inequality 

(l - 1) < (l + I) 5" + ^ll^ll^. (96) 

We assume that 

T 

1 = exp(— 0.5t), 

and therefore 

1 + - < exp(0.5r). 

£ 

With our notation it is easy to obtain from (|96|) the required estimate of stability 

(EH). 

This estimate of stability is characterized by using more complex norms in 
compare with the case of standard schemes with weights (compare (l35t and 
(|9"Tj) ). The numerical implementation of the regularized scheme (I5i?|) . (j^H)) can 
be performed similarly to scheme (I75p . (|76[) . 

8. Numerical results for model problems 

Numerical experiments are performed for the parabolic equation where 

fc(x) = 1, /(x, t)=0, x 6 O, < i < T. (97) 

Problem (f9|)— (fTT|) . (|96| is considered in the unit square ^ = / 2 = 1, and the 
initial condition has the form 

w°(x) = sin(ni7rxi) sin(ri27TX2), x € f2, (98) 

for natural rti and n-i. The solution of problem (|9j)- (fTTj) . ((97]), ((98]) is written 

as 

u(x, t) — exp(— ir 2 (n 2 + n\)t) sin(ni7ra;i) sin(n27ra2). (99) 
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Figure 5: The error of factorized regionally-additive scheme 



The numerical results derived using the regionally-additive schemes are com- 
pared with the difference solution obtained by means of the implicit scheme (|30|) , 
(p{Tj) for a = 1/2 and a = 1. The error of the approximate solution was estimated 
as e(t n ) — ||y"(x) — u(x, t n )\\ at each particular time level. 

In the basic case we used ni — 2,n 2 = 1, ATj = iV 2 , hi = h 2 = h = 1/40, 
T = 0.05, N = 10, t = 0.01. The decomposition is carried out by cutting the fi 
into four squares (h — 0.5). 

Results of solving the test problem using the standard implicit schemes with 
weights ([30]) . ([3i~]) of second {a = 1/2) and first (a = 1) orders of accuracy with 
respect to r and factorized regionally-additive scheme (|39]l. ([54]) . (|55|) (FAS) for 
the same values of the weight parameter a are presented in FigJS] 

With the selected parameters the domain decomposition scheme, constructed 
on the basis of classical factorized schemes, yields the approximate solution 
with a slightly larger error than the standard two-level scheme with weights. 
The effect of time step is presented in Fig. [51 where data are obtained with 
redusing time step (r = 0.005). A more interesting effect is connected with the 
discretization in space (Fig|7]). The effect of conditional convergence becomes 
more evident for the regionally-additive scheme at a = 1/2. 

In the study of the decomposition schemes particular attention should be 
paid to the dependence of accuracy of the approximate solutions on the number 
of subdomains. The error of the schemes for the increased number of subdomains 
(four times (h = 0.25)) is shown in FigJS]). Decreasing of the accuracy (compare 
Figj5]and Figj8]) is more significant for the factorized regionally-additive schemes 
at g — 1/2. 

For the model problem ©-CH]), (J9TJ) , ((981) witri tne exact solution <(99j) and 
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Figure 6: Reducing of the time step (t = 0.005) 



decomposition (|40)) . (|45|) we performed calculatuins via the above schemes of 
multicomponent-wise splitting. The error of the approximate solution for the 
basic case derived using the regularized regionally-additive scheme (|75|) . (|76j) 
with cr = 1 is compared with the results of regionally-additive component- 
wise splitting scheme (|55|) with a = 1/2 in FigO The accuracy of the 
regularized scheme is clearly much lower. As for the accuracy of the component- 
wise splitting scheme, it is practically the same as the accuracy of the factorized 
regionally-additive scheme (see FigEJand FiglBj). 
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Figure 7: Reducing of the spatial step (h = 1/80) 
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Figure 8: Increasing of the number of subdomains (h = 0.25) 
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